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Abstract 

We derive exact results for gap probabilities, as well as densities of extreme eigenvalues for six complex random 
matrix ensembles of fundamental importance. These are Gauss-Wigner, Laguerre-Wishart, Cauchy-Lorentz 
(two variants), Jacobi-MANOVA and Bures-Hall (trace unrestricted). We deal with both correlated and uncor¬ 
related cases. Extensive Monte-Carlo simulations based on explicit matrix models or Dyson’s log-gas formalism 
have also been performed, which corroborate all analytical results. 
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1. Introduction 


Given a random matrix ensemble, the question regarding probabilities of gaps in its eigenvalue spectrum 
is of fundamental importance. The investigation of gap probabilities is not only interesting from a general 
mathematical point of view, but also because of its connection with several interesting problems and concrete 
applications. Naturally there has been a plethora of studies which deal directly or indirectly with this ques¬ 
tion [i- 


. For example, in the case of a Hessian matrix drawn from a Gaussian ensemble, which can model a 


random multi-field potential or a landscape, the gap probability serves as an estimator for the fraction of local 
minima amongst the stationary points of the landscape [ll-y|. As application to geometry and random topol¬ 
ogy, gap probability has been used to calculate the intrinsic volume of the set of Frobenius norm one singular 
random matrices, and computation of certain class of Betti numbers P|. In Q gap probability has been shown 
to be related to the asymptotics of orthogonal polynomials defined on an arc of unit circle. In Q| a connection 
has been established with the expected number of minima of a random polynomial. Gap probabilities are also 
relevant in deciding physical stability in dynamical systems and ecosystems HB- 

Intimately related to gap probabilities is the topic of extreme eigenvalues, which has been of considerable 
interest to researchers in several fields [3|, llSHSSj . One of the notable examples is the celebrated Tracy-Widom 



distribution for extreme eigenvalues mm I, which emerges in seemingly unrelated problems, such as mesoscopic 
fluctuations of excitation gap in quantum dots 0 10 , distribution of the pseudo-crtical temperature in mean- 


field spin glasses [421 


sequence alignment 


height fluctuations of non-intersecting Brownian motions 


m 


44l | , growth models 45|-|48| , 


491. random permutations [0 . finance M|, etc. In the estimation of performance of 
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multiple channel communication systems the statistics of extreme singular values (square root of eigenvalues) 


of channel matrix plays a crucial role 


19-3- In the problem of quantum entan glem ent the extreme Schmidt 


eigenvalues carry important information about the degree of entanglement 
portfolio is related to the extreme eigenvalues of the correlation matrix 
determination of metric such as condition number 


53 


24 


26l |. In finance, the optimal 


52l| . Extreme eigenvalues also lead to 


54| . which is a measure of the relative conditioning 


(or rank-deficiency) of a matrix and has numerous applications in a variety of contexts, for example multiple 
channel communication systems, linear detection and classical linear algebra. Moreover, in r ecent y ears the large 
deviations of extreme eigenvalues of random matrices have been studied in several contexts 


55 


59l| which range 

from principal component analysis to Wigner time-delay distribution in chaotic cavities. A recent survey can 

62|. 


be found in [3 ■ Large deviation studies relate to yet another riveting topic of extreme or rare events 


61 


The case of matrices with large dimensions has been exhaustively studied for extreme eigenvalues, while for 
finite dimensionality case there have been relatively fewer studies with analysis remaining concentrated mostly 
around Laguerre-Wishart ensembles. Our aim here is to derive finite dimensionality exact results for gap 
probabilities, and probability density functions for extreme eigenvalues. We cover the following ensembles of 
complex matrices which are of central importance to random matrix theory and have wide range of applications: 


1. Gauss-Wigner Ensemble 

2. Laguerre-Wishart Ensemble 

3. Cauchy-Lorentz Ensemble (Variant I) 

4. Cauchy-Lorentz Ensemble (Variant II) 

5. Jacobi-MANOVA Ensemble 

6. Bures-Hall Ensemble (trace unrestricted) 


We consider both correlated, and uncorrelated cases of these ensembles. In the case of uncorrelated ensem¬ 
bles the underlying matrix probability measure is rotationally invariant. On the other hand, associated with 
correlated ensembles are rotationally noninvariant probability measures which are comparatively difficult to 
handle. As we find below, these ensembles fall broadly under two categories: the ones with joint density of 
eigenvalues involving product of two determinants, and the others with joint density involving product of a de¬ 


formula 


64 


terminant and a Pfafhan 3 • Therefore, we also provide general results based on generalizations of Andreief’s 


66 1 and de Brujin’s formula 671. l68|. 


The remaining of the paper is organized as follows. We focus on the generic setup in Sec [2] and derive 
general results for the two categories of ensembles indicated above. Sections [3] to |8] deal with explicit results 
for the six (times two) random matrix ensembles listed above. We conclude in section 12] with a brief summary 
and outlook. 
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2. General results 


2.1. Preliminaries 

Consider the joint probability density function (JPDF) P({A}) of n eigenvalues associated with 

certain random matrix ensemble, such that 


/ dAi-- - / dA„P({A}) = l. 

JO JO 


( 1 ) 


Here II and tu are the lower and upper integration limits and set the domain for the eigenvalues. In the cases 
when P({A}) has an associated easy to implement explicit matrix model, one can generate the eigenvalues 
by diagonalization. While for any well behaved JPDF, with or without any explicit matrix model, one can 


always generate the eigenvalues using the Dyson’s log gas formalism 


69 


70|. Given this, one is naturally 


interested in some metric to determine the fraction of eigenvalues which lie in a certain interval within the 
domain of the eigenvalues. This is equivalent to figuring out the probability that a certain interval is devoid 
of any eigenvalues, and hence one is led to the study of gap probability. A closely related problem is the 
investigation of behavior of the smallest or the largest eigenvalue, for which the full information is provided 
by the corresponding probability densities. In this section we formally define these quantities and afterwards 
work out the expressions for two kinds of structures for the matrix ensembles which are very generic and cover 
a number of important cases. 

We will calculate the gap probability 


A(r, s) = J 


( 2 ) 


dXi--- J dA„P({A}), 

(0,r)U(s,iu) (D,r)U(s,fu) 

which refers to the probability of finding no eigenvalue in (r, s), such that II < r < s < (u. We also examine the 
probability that all the eigenvalues lie in [r, s]: 


E{r, s)= dXi 


dXnP{{X}). 


( 3 ) 


We should underline that 


E{r,s) ^1-E{r,s), (4) 

except when n = 1. We observe that E{r,s) is same as the gap probability that there are no eigenvalues in 
either of the intervals (0, r) and (s, u). We will refer to it as the double gap probability. 

Of particular interest are the gap probabilities E{i,x) and E{x,u), with 0 < cc < (u. Clearly 


E{i,x) = E{x,u) = / dXi--- / dXnP{{X}) 


( 5 ) 


gives the probability that there are no eigenvalues between 0 and x, which is same as the probability that all 
eigenvalues are greater than or equal to x. As such, it refers to the survival function (SF) or reliability function 
of the smallest eigenvalue. On the other hand 


E{x,u) = E{i, x) = / dAi • • • / (iA„P({A}) 


( 6 ) 
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gives the probability that there are no eigenvalue between x and u or that all eigenvalues are less than or equal 
to X. Therefore, it is the cumulative distribution function (CDF) of the largest eigenvalue. Survival function 
and cumulative distribution function are related as SF= 1—CDF. 

The probability density function (PDF) of the smallest eigenvalue, and that of the largest eigenvalue are 
obtained respectively as 


and 


Ps{x) = - 


Pl(x) = 


^^^(O,^) dE{x,<ii) 

dx dx 

dE{x,<ii) dE{},x) 


( 7 ) 


( 8 ) 


dx dx 

We may also evaluate the joint probability density of the smallest and the largest eigenvalues using E{r, s) as 


Psl(7’, s) = - 


d^E{r, s) 
drds 


( 9 ) 


With the information of psl(?') s) "we can examine the statistics of quantities which involve the smallest eigen¬ 
value as well as the largest eigenvalue. For instance, we can compute the density of condition number. 

As already indicated, in the matrix ensembles that we pursue in the following sections, the JPDF exhibits 
two kinds of structure. The first one has a biorthogonal form of Borodin type |7l| and involves product of two 
determinants. The second one involves, on the other hand, the product of a determinant and a Pfafhan. We 
will refer to these as ensembles of Type I and Type II, respectively. We present the general results for these 
ensembles below, and explicit results for the particular cases are given in the subsequent sections. 


2.2. Type I Ensembles 

We consider in this category, the ensembles which involve product of two determinants, namely 




( 10 ) 


1=1 


Here we use the notation | ■ | to represent determinant. We note that the product involving the weight function 
w{X) can be absorbed within the determinants. This is clearly of Borodin’s biorthogonal form [71|. We will 
deal with correlated ensembles for which we find fj{Xk) = which makes the corresponding determinant 
the Vandermonde determinant: 


A„({A}) = |Ar'|,.fc=i.....n = 


( 11 ) 


j>k 


Eq. (HOI) also includes the case of classical unitary ensembles which arise when in addition to /j(Afc), gj{Xk) = 
as well. This situation is encountered in the case of uncorrelated ensembles. 


Ar' 


For this class of ensembles Andreief’s integration formula 
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66| yields at once the expression for the 


normalization factor, as well as the gap probabilities E(r,s) and E(r,.s). The inverse of normalization factor 
(partition function) is found to be 

= n\ \hj^k\j,k=i,...,n (12) 
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where 


hj^k= / dXw{X)fj{X)gkiX) 


(13) 


The gap probability is obtained as 


E{r, s) = nl C\xj,kir, s)\j,k=i,...,r 


(14) 


with the kernel Xj,/c(r, s) given by 

Xj,k{r,s)= J dXw{X)fj{X)gkiX) (15) 

(D,r)U(s,aj) 

The expressions for i?(ll, s) and E{r, u) can be read easily from Eq. (jl4p with the kernel Xj,k getting simplified 
respectively to 

Xi,fc(M) = / dXw{X)fj{X)gkiX), (16) 

J S 

and 

Xj,k{r,'»)= j dXw{X)fj{X)gkiX). (17) 

Similarly, E{r,s) has the expression 


E{r, s) = nl C'|xj,fc(r, s)\j^k=i,. 


(18) 


where the kernel Xj,k{f, s) is given by 


Xj./c(r,s)=/ dXw{X)fj{X)gk{X). 
J r 


(19) 


The expression for density of the smallest eigenvalue is obtained using Eq. ([7]) upon differentiating the 
determinantal expression as 


ps{x) =n\CY^ \(l)fl{x)\j,k=i 


( 20 ) 


where 


(i) . , J j = i, 


= 


( 21 ) 


Xj,ki^,x), j^i. 

In a similar manner Eq. (|8]) gives the density of the largest eigenvalue as 


Pl{x) =n\Cj2 \'^jlix)\j,k=i,..., 


( 22 ) 


with 


X _ J ^(2;)/j(a;)5fe(a;), j = i, 


= 


(23) 


^Xj.fe(a;.'u), j¥^i- 

It is to be noted here that Xo,k{^}^) = Xj,k{x,u), and Xi,k{x,u) = Xj,ki^Tx). 
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2.3. Type II Ensembles 

We consider in this category the ensembles whose joint probability density involve product of a Pfaffian 
(Pf) and a determinant, i.e. 




(24) 


1=1 


The kernel /(A,/i) is antisymmetric in A and i.e., = —f{p,X). Ensembles of above type have 

been encountered in the study of Wilson random matrix theory . Moreover, the JPDF of eigenvalues for 
unrestricted trace Bures-Hall ensemble can also be cast in the above form [7^; see section |8] ahead. The 
normalization factor, as well as the gap probability expressions can be obtained using a generalization of de 
Bruijn integration theorem 
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63. 


Let us define for convenience 

i, n even, 
n + 1, n odd. 
We then have the expression for the partition function as 


N = 




(25) 


(26) 


where 


hj,k= d\ dpw{X)w{p)f{\,fj.)gj{X)gk{fJ.) 
Jw Jw 


when n is even, and in addition 


dX / dp,w{X)w{p.)f{X,p) [gj{X)gk{g) - gk{X)gj{p,)] 


^dJ 

hj^n+i = -hn+i,j = (1 - Sj^n+i) / dXw{X)g{X) 

Jn 


(27) 


(28) 


when n is odd. Here, Sij is the Kronecker delta and therefore the factor (1 — (5j_n,_|_i) ensures that hn+i^n+i = 0. 
We found that the second expression in Eq. (l27t is more stable if used for numerical evaluation. 

The gap probability can also be written in terms of a Pfaffian as 


E(r, s) = n'.CFi[xj,kir,s)]j^k=i,...,N, 


(29) 


with the kernel Xj,k{r, s) given by 


Xj.k{r,s)= J dX J 


dp. w(X)w{p)f{X, p)gj {X)gk{p) 


(0,r‘)U(s,(u) (ll,r)U(s,dj) 

1 . 


I dX J dpw{X)w{p)f{X,p)[gj{X)gk{p) - gkiX)gj{p)] 

(D,r)U(s,aj) (D,r)U(s,iu) 

when n is even. When n is odd, additionally we have 

Xj,n+iir,s) = -Xn+ij{r,s) = (1 - Sj^n+i) J dXw{X)gj{X). 

(0,r)U(s,iu) 


(30) 


(31) 
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Similarly, the double gap probability also turns out to be a Pfafhan: 


E{r,s) = n\CPt[xj,k{r,s)]j^k=i,...,N- 

Here the kernel s) is 

Xj,k{r,s)= dX dfj.w{X)wip)f{X,fi)gj{X)gkitJ-) 

Jr Jr 


r Jr 

s pS 


dX / dg.w{X)w{n)f{X,fi)[gj{X)gk{g.) - gk{X)gjig.)] 


/ r Jr 


when n is even, and in addition 


Xj,„+i(r,s) =-x„+ij(r,s) = (1 - (5y„+i) / dXw{X)gj{X) 

J r 


(32) 

(33) 

(34) 

(35) 


when n is odd. 

The expression for density of the smallest eigenvalue is obtained using Eq. dZD by differentiating the Pfaffian- 
expression as 


1 ^ 

Ps(.x) = -n\c Pi~^[xj,k{Kx)]j,k=i,-,N^\ 4 >%{x)\jx=i,...,N, 


with 


Likewise, the expression for density of the largest eigenvalue is given by 

1 

pjx) = -n\c Pi~^[xjx{x,'u)]jx=i,...,N^\i^fu^x)\jx=i,...,N, 


(36) 


(37) 




where 


Xj,k{x,u), j^i. 


(38) 


(39) 


We now specialize to specific random matrix ensembles and work out the explicit results using the above 
general expressions. 


3. Gauss-Wigner ensemble 

3.1. Correlated case 

We consider n-dimensional Hermitian matrices H from the distribution 

T’(H) oc (40) 

where, without loss of any generality, we take S = diag (cti, ..., (T„), with all cr’s greater than zero. In case S 
is not already in diagonal form, tr’s are the corresponding eigenvalues. Similar consideration applies to other 
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Table 1: Gap probabilities: Comparison between analytical and simulation results for correlated Gauss-Wigner 
ensemble. In each case the a values are from (tri, cr 2 , 0 - 4 , cts) = (3/4,4/9,1/4, 2/7, 5/6). 


n 

r 

s 

E{r,s) 

E{r,s) 

Analytical 

Simulation 

Analytical 

Simulation 

2 

- 0.5 

0.6 

0.0808 

0.0804 

0.3286 

0.3281 

3 

0.1 

1.4 

0.0347 

0.0346 

0.0017 

0.0018 

4 

— CX) 

-0.2 

0.0111 

0.0112 

0.0000 

0.0000 

4 

-0.38 

0.38 

0.0034 

0.0031 

0.0012 

0.0011 

5 

-0.7 

00 

0.0000 

0.0000 

0.1541 

0.1537 


correlated ensembles solved ahead. From Eq. (1401) it is evident that the matrix elements of H are zero-mean 
Gaussians having variances as 

^2 


(Hf,) = 


((ReH,fe)') = ((ImH.fc)") = 




(diagonal), 

, (off-diagonal). 


(41) 


2(^1+ 

We refer to this case as a correlated Gauss-Wigner ensemble. 

The joint eigenvalue density for Eq. (Hill) can be obtained with the aid of Harish-Chandra-Itzykson-Zuber 
(HCIZ)-integral [t^i 3 written down as 


^({^}) “‘^A+({A})I kfc=i. -.- 


(42) 


Here A({A}) = ~ ^k) is the Vandermonde determinant, and A-|_({A}) = + ^k)- The 

eigenvalues lie in (— 00 , 00 ). Schur’s Pfafhan identity gives 72l.I75ll76l|. 


- ^k ^ 


Pf [(Afe - Xj)/{Xk + Aj)]. 




Pf 


[(Afc [l]i=l,...,n 


n even. 


odd. 


(43) 


As we can see, this ensemble falls in type II according to our classification, with w{X) = 1, fj{Xk) = e 
and g{Xj,Xk) as in Eq. (l43l) . The normalization factor in Eq. (l42l) can be obtained using Eq. (IMl) as 


n! Pf [nakajial - a^)/{(Tl + , 


c-^ = 


i!Pf 


[nakajial - af)/{al + 


[ V^'^k]k=l 


n even, 


odd. 


(44) 


= n) 


i—l j>k ^ 




(45) 



































X X 

Figure 1: Plots for correlated Gauss-Wigner ensemble with n = 3, and (cti, cr 2 , ca) = (3/5, 2/3, 5/4). (a) SF of 
the smallest eigenvalue (b) PDF of the smallest eigenvalue (c) CDF of the largest eigenvalue (d) PDF of the 
largest eigenvalue. 


The gap probability is given by Eq. (1^ with 


and 


XjMr,s)= j 


dX 




( —oo,'r)U(s,oo) ( —cxD,r)U(s,oo) 

1 

■ 2 


/i + A 


e ■> e 


I I 


Li — A / —2^2 —2 2 -2 2 —2 ^2 

P / —A —G, a —G- a —G, A 

‘ e J e ^ ^ — e ^ ^ e ^ 


(46) 


( —oo,r)U(s,oo) ( —oo,r)U(s,oo) 


y — ^ \ ^ 

dX e 

( —oo,r)U(s,oo) 

= [2 + erf ((T“ V) - erf {aj^s)] (1 - 5j,„+i), 


(47) 


(48) 


where erf( 2 :) = {2/x/n) fgdte is the error-function. Similarly Xi.fe(r, s) gives the double gap probability via 
Eq. (15^ with 
































































and 


Xj,n+i{r, s) = -Xn+1,j{r, s) = (1 - Sj^n+i) / dXe 


= —^3 [erf (ctj ^s) - erf {a^ V)] (1 - 5^,„+i). (50) 

The expressions for density of the extreme eigenvalues follow from Eqs. dSH) and (1551) . 

In Table [T] we compare the results for gap probabilities for various choice of parameters. The simulation 
results have been calculated using 50000 matrix realizations, as has been also done for other ensembles ahead. 
Fig. m shows the plots of the distribution functions and densities obtained from analytical expressions (solid 
lines) as well as from Monte-Carlo simulations (histograms). We find excellent agreement in all cases. We 
note that there is a symmetry x(—oo,—a;) = x(a;, oo) in this case. It leads to E{—oo,—x) = E{x,oo) and 
ps{—x) = pl{x) which are clearly observed in the plots. These correspondences are consequence of the 
symmetry exhibited by the matrix model (1401) . 


3.2. Uncorrelated case 

We consider here the case S = !„ in (l40l) . and hence obtain 


iP(H) oc 


which is the density for the classical Gaussian Unitary Ensemble (GUE) 


(51) 
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771 . In this case the variance of 


diagonal elements of H is 1/2, while that of the real and imaginary parts of off-diagonal elements is 1/4. The 
results for the present case can be obtained from the preceding section using a limiting procedure to set all a 
equal to 1. However, it is more straightforward to start with Eq. m, which gives the following result for the 
joint density of eigenvalues: 

n 

P(Ai,..., Xn) = GA2({A}) n (52) 

3 = 1 

It is interesting to see how the product of a Pfaflian and determinant gives rise to the product of two 
determinants for S — In. Gomparing Eq. dsa) with our type I ensemble form, Eq. we find that 

fj{Xk) = Qi(X k) = Xir^, and weight function w{X) = . The normalization factor in Eq. (1551) is well 

known [631 . l77|| : 


C-^ = nl 


i + (-iy** p/j+t-i' 


,(^) 


n/2 




nro+i)= 


7r"/2 


3 = 1 


2n(n-l)/2 


G(ti + 2), 


(53) 


where r(z) = dtU~^e~* is the Gamma function, and G{n) is the Barnes G-function which is defined for a 
positive integer n as G(n -b 1) = 0^=1 r(/). 

Gaussian ensembles are probably the most studied random matrix ensembles, along with the Laguerre- 
Wishart ones 63], l77|. For large n, the distributions of extreme ei genv alues for Gaussian ensembles, including 


the complex case above, were derived by Tracy and Widom in 


13 


14| . Further progress has been made 
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Table 2: Gap probabilities: Comparison between analytical and simulation results for uncorrelated Gauss- 
Wigner ensemble. 


n 

r 

S 

Eir, s) 

E{r,s) 

Analytical 

Simulation 

Analytical 

Simulation 

2 

-2 

-0.2 

0.1814 

0.1811 

0.0418 

0.0421 

3 

0.05 

CX5 

0.0074 

0.0075 

0.0042 

0.0041 

3 

— CX) 

0 

0.0056 

0.0056 

0.0056 

0.0053 

4 

-1.2 

1.2 

0.0018 

0.0019 

0.0073 

0.0073 

5 

-1.4 

4 

0.0000 

0.0000 

0.0589 

0.0587 


in [171, l55|, l58| where large deviations of the extreme eigenvalues have been explored in detail. We are interested 
here in finite n exact results. 

The expression for gap probability is given by Eq. m with 


(1 + (_iy+fc) r(i±|^) + (-iy+'= r2) - s")) 


X3,kir,s)= i (_l)J+fcr(i±|^,r2) + r(i±|^,s2) 


1 

2 


r < s < 0, 
r < 0 < s, (54) 
0 < r < s. 


(1 + (_i)j+fe)r(2±|^) - r(i±|^,r2) + r(i±|^,s2) 

Note that when j' + fc is odd, all the three expressions above become identical. The expression for double gap 
probability E{r, s) is given by Eq. (1181) with 

(3+k-l „2 


XjAes) = < 


i(_i)j+fc _r(2±|^,r2) + r(^±|^,s2) 


(1 + (_i)j+fe) - (-i)j+fcr(i±|^, r2) - s2) 


r(i±|^,r2)-r(^,s2) 

These reduce to the following for r = —oo or s = oo: 

'(1 + (-i)J'+'=)r(i±|^) - (-i)j+'=r(i±|^,x2) 


Xj,ki-oo,x) = Xj,k{x,oo) = 


ir(i±^,x2), 


1 (_i)j+fer(4±|^,x2) 


r < s < 0, 
r < 0 < s, 
0 < r < s. 


X < 0, 
X > 0, 


Xj,fe(x,oo) = Xj,fe(-oo,x) = 

The density of smallest eigenvalue is given by Eq. dini) with 

= 

and that of the largest eigenvalue by Eq. (1^^ with 


(1 + (_i)j+fc) r(i±|^) - r(i±^, x2) 


X < 0, 
X > 0. 


e ^ j = i, 


Xi,fc(-oo,x), j^i. 


V'Jka;) = 


e j = i, 


Xi,fc(a;,oo), j^i. 


(55) 


(56) 


(57) 


(58) 


(59) 
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Figure 2: Plots for uncorrelated Gauss-Wigner ensemble with n = 6. (a) SF of the smallest eigenvalue (b) 
PDF of the smallest eigenvalue (c) CDF of the largest eigenvalue (d) PDF of the largest eigenvalue. 

In Table [2] we compare the results for gap probabilities for various combination of parameters. Fig. [2] 
shows the analytical distributions and densities, as well as histograms obtained from Monte-Carlo simulations. 
We again find the symmetries E(—oo,—x) = E(x,oo) and ps{—x) = pl{x) manifest in the plots. We note, 
however, in this case that Xj,k(,—oo, —x) = {—iy~^^Xj,k{x, oo) and the aforementioned symmetries follow from 
the determinanatal expression m- The underlying reason for these is again the matrix probability density. 


4. Laguerre-Wishart ensemble 

4.I. Correlated case 

The correlated Laguerre-Wishart ensemble is described by the density 

iP(H) oc (60) 

for positive-definite (or more generally nonnegative-definite) Hermitian matrices H with a > — 1 and S = 
diag (cti, ..., cr„), tTj > 0. For nonnegative integer values of a this ensemble can be realized as follows. Consider 
n X m-dimensional (n < m) complex Gaussian matrices G from the distribution 

iPG(G) oc (61) 

then n-dimensional positive-definite Hermitian matrices (Wishart matrices) 

H = GG^' (62) 
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constitute the correlated Laguerre-Wishart ensemble with the associated density given by Eq. dSO]), and the 
parameter a = m — n. We underline that the m x m matrices G^G have m — n generic zero eigenvalues, while 
the n nonzero eigenvalues are same as those of GG^. The parameter m is referred to as the degree of freedom 
of the Wishart distribution. Equivalent to the above construction, if n x m-dimensional (n < m) complex 
(Ginibre) matrices G are taken from 

iPg(G) oce-*’-®®', (63) 

then H of Eq. (l60ll with a = m — n can be generated as 

= (64) 


The joint eigenvalue density for the correlated Laguerre-Wishart ensemble can be obtained with the aid of 
HCIZ-integral, and written in a biorthogonal form as 


P{Xu-.,Xn) = CA{{X})l[Xf 


I ^ — ^ At. I 


(65) 


i=i 


with 0 < Aj < oo. This JPDE falls in type I as per our classification with w{X) = A“, fj{Xk) = Xj. ^ and 
9j{Xk) = e~'^i . The partition function for Eq. (1^ is 


C-^ = n\ r(j + , 




= n\A{{a})l[af+^T{j + a). 


( 66 ) 


i=i 


Eq. (I62|) or (|Mll provides a way to easily generate matrices and hence the eigenvalues distributed according to 
Eq. (|65l) for non-negative integer a values, while for any real a (> — 1) the eigenvalues can be generated with 
the aid of Monte-Carlo simulation based on Dyson’s log-gas picture [69|, l70| , as already indicated in Section [5J 
Laguerre-Wishart ensemble has been explored extensively because of their crucial role in the field of multi¬ 


variate statistics 


telecommunication 


79l|. a nd explicit appearance in problems related to time series 


m 1 J 23 


SOj-lS^ and multiple-channel 


83l-l85l|. As a consequence several results concerning extreme eigenvalues are avail- 


34 


able. Some of the most recent results being due to Wirtz el al. 
PDE of extreme eigenvalues possess form as in 
and closed structure. 

The gap probability in this case is given by Eq. (fT4l) with 


Sj-^. Our results for gap probabilities and 


35l| . however the kernels involved have relatively simpler 


XjAas) = + a,a^ V) -h r(j -b a,a^ ^s)], ( 67 ) 

where 7 ( 0 , r) = dz and r(a,r) = dz z°'~^e~^ represent the lower-incomplete gamma function 

and the upper-incomplete gamma function, respectively. Also, for the double gap probability, as in Eq. (fT51) . 
we have 

Xi.fc(As) = aA°‘[T{j + a,aAr)-T{j + a,aAs)]. (68) 

These simplify to the following for r = 0 or s = c»: 


Xi.fc( 0 , x) = x{x, 00 ) = cr^’^“r(j -b a, a^. ^a;). 


(69) 
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Table 3: Gap probabilities: Comparison between analytical and simulation results for correlated Laguerre- 
Wishart ensemble. The a values are from (cti, cts) = (2,3/5, 7/3, 2/5,4/9,1/2). The m values are indicated 
when matrix construction is possible using Eq. (EH) or 


n 

a 

m 

r 

s 

E{r, s) 

E{r,s) 

Analytical 

Simulation 

Analytical 

Simulation 

2 

1.25 

- 

0.63 

5 

0.1789 

0.1786 

0.1913 

0.1926 

2 

3 

5 

3.1 

CX) 

0.0046 

0.0047 

0.1753 

0.1747 

3 

2 

5 

0 

5 

0.0013 

0.0013 

0.0006 

0.0007 

4 

4 

8 

2.5 

20 

0.0020 

0.0020 

0.0103 

0.0105 

5 

-0.3 

- 

4 

24 

0.0051 

0.0053 

0.0000 

0.0000 

6 

0 

6 

0.3 

28 

0.0000 

0.0000 

0.0523 

0.0521 



Figure 3: Plots for correlated Laguerre-Wishart ensemble with n = 3 ,to = 7, and S = (1/2, 3/4,4/5). (a) SF 
of the smallest eigenvalue (b) PDF of the smallest eigenvalue (c) CDF of the largest eigenvalue (d) PDF of the 
largest eigenvalue. 


Xyfc(a;,oo) =x{0,x) =cr^+“7(j + a,crj,^a;). 


(70) 


The expressions for the smallest eigenvalue and the largest eigenvalue are given respectively by Eqs. (1201) 
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and (|2^ with 


and 




V'S(a;) = 



j = i. 

x(0,a;), 

j ^ *■ 

2.J+a-l 

j = i 

x{x,oo), 

j 7^ i 


(71) 


(72) 


In Table |3] we list the gap probabilities for several choices of parameters. Fig. |3] displays the plots of dis¬ 
tribution functions and probability density functions obtained from analytical expressions, which are validated 
by the histograms obtained from Monte-Carlo simulation. 

A particularly interesting scenario occurs when a = 0, i.e., m = n. In this case the gap probability E(0,x) 
or E{x,oo) and the density of the smallest eigenvalue possess remarkably simple expression. These can be 
obtained conveniently if we start from the JPDF in Eq. (I60|) and then implement the translation = Xj — x. 
Upon a little simplification, the gap probability is obtained as 


E{0, x) = E{x, oo) = exp I —x ^ ^ ■ 

V j=i 

The density of the smallest eigenvalue then follows immediately as 


(73) 


ps{x) = M 




(74) 


u'=i 


1=1 


4-2. Uncorrelated case 

We now deal with the uncorrelated Laguerre-Wishart case which involves positive-definite Hermitian ma¬ 
trices H from the density 


P(H) oc |H|' 


„-trH 


(75) 


with a > — I. Again, for nonnegative integer values of a (= m — n) the matrices from the above density can be 
generated as H = GGl, where G are nxm dimensional (n < m) complex matrices taken from the density (1631) . 

Equation (123 is recognized as the standard Laguerre Unitary Ensemble of random matrices with the 
corresponding joint eigenvalue density |63l. l77l| : 


P(Ai,...,Aw) = CA2({A})nA° 




(76) 


i=i 


This is ensemble of type I according to our classification. When compared with Eg . (HQl) we find that fjiXk) = 
gj{Xk) = Xi~^, and w(X) = A“e“^. The normalization factor is well known 63j, [TTf: 


C = n<.\r(j + k + a - 

n 

= nr(j + i)r(j + a)- 

i=i 


(77) 
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Table 4: Gap probabilities: Comparison between analytical and simulation results for uncorrelated Laguerre- 
Wishart ensemble. The m values are indicated in the cases when matrix construction is possible using GG^ 
where G is from Eq. (1551) . 


n 

a 

m 

r 

s 

E{r, s) 

E{r,s) 

Analytical 

Simulation 

Analytical 

Simulation 

2 

2 

4 

1 

5 

0.1221 

0.1224 

0.2298 

0.2300 

3 

-0.25 

- 

0.9 

9 

0.0085 

0.0088 

0.0273 

0.0269 

4 

0 

4 

1.5 

14 

0.0001 

0.0001 

0.0021 

0.0020 

5 

2.6 

- 

0 

15 

0.0000 

0.0000 

0.2554 

0.2546 

6 

3 

9 

18 

CX) 

0.1643 

0.1640 

0.0000 

0.0000 


The expression for gap probability is given by Eq. (j29p with the kernel 

s) = 7(j + fc + a - l,r) + r(j + fc + a - l,s). 

Similarly, the double gap probability is given by Eq. (15^ with 

Xi,fc(G s)=T{j + k + a- 1, r)-T{j + k + a- 1, s). 

For r = 0 or s = 00 these yield 

Xj,fc(0,x) = Xj,k{.x,oo) = T{j + k + a- l,a:), 

Xj.fc(a:,oo) = Xj,fc(0,x) = 7(j + fc + a - 1, x). 

The expressions for densities of the extreme eigenvalues are given by Eqs. (1201) and (l22l) with 

aO) 


and 


= 




j.j+k+a-2 g-x 

II 

X(0,x), 

j + 

^j+k+a-2 g—X 

j = * 

x(a;,oo). 

j ^ i 


(78) 

(79) 

(80) 
(81) 

(82) 

(83) 


Table [4] displays the results for gap probabilities with comparison between analytical result and numerical 
simulation. Fig. 2] shows the plots for survival function for the smallest eigenvalue, cumulative distribution 
function for the largest eigenvalue, and the corresponding densities. 

Following the correlated case, Eqs. dZSl) and dill), for a = 0 (to = n) the gap probability E(0, x) or E{x, oo) 
and the density of the smallest eigenvalue have the following simple expressions [l^ : 

x) = E{x, oo) = (84) 

p$(x) = ne“”“. (85) 
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Figure 4: Plots for uncorrelated Laguerre-Wishart ensemble with n = = b. (a) SF of the smallest eigenvalue 

(b) PDF of the smallest eigenvalue (c) CDF of the largest eigenvalue (d) PDF of the largest eigenvalue. 


5. Cauchy-Lorentz ensemble: Variant I 

5.1. Correlated case 

We refer to the ensemble of n-dimensional Hermitian matrices H from the density 

iP(H)oc |1„ + S-2 h2|-", (86) 


as the correlated Cauchy-Lorentz (variant I) ensemble. Here S = diag (cti, ..., cr„), aj > 0, and k > n — 1/2 for 
convergence. Using unitary group integral result of we arrive at the joint eigenvalue density 


PiW) = 




(87) 


A+({A})' 

with —00 < Xj < oo. Similar to the correlated Gauss-Wigner case, the above JPDF can be rewritten as product 
of a determinant and a Pfaffian. Therefore, the normalization factor is provided by a Pfaffian expression as 


C = Fi[hjk]j,k=i,.., 

where N is as defined in Eq. (1251) . Also, for even n we have 
7rr^(K — n + 1/2) 


N, 


hj,k — 


2T^{K-n + l) 


^ 2 (k n+i)^^2(K ") ^2 (k — n)-I-1, K — n-b 3/2; 2 (k — n-b 1); 1 — tTfe 

(2{k — n) + 1, K — n + 3/2;2 {k — n + 1);1 — . (88) 


2(k —n+1) —2 (k— n) 7-, 

-ctU V, '2F1 
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Table 5: Gap probabilities: Comparison between analytical and simulation results for correlated Cauchy- 
Lorentz (variant I) ensemble. The a values are from (cti, erg) = (5/6, 7/8,16/11, 2,3/10). 


n 

K 

r 

5 

Fir, s) 

E{r, s) 

Analytical 

Simulation 

Analytical 

Simulation 

2 

3 

-3 

0 

0.1251 

0.1242 

0.1078 

0.1073 

3 

3.2 

-0.5 

1.2 

0.0471 

0.0492 

0.0116 

0.0118 

3 

4 

-0.6 

0.7 

0.0356 

0.0354 

0.0195 

0.0198 

4 

4.5 

2.6 

16 

0.6117 

0.6108 

0.0000 

0.0000 

5 

5 

-9.5 

9.5 

0.0000 

0.0000 

0.6423 

0.6440 


and when n is odd, additionally we have 

/ OO 

dA (1 + A ) 

-OO 

r(«; — n + 1/2) 


— 2 \ 2'^—K--t-n—1 


= y/lTO 


■(1 Sj^n+l)- 


(89) 


^ r(K-n + l) 

In the above expression for hj^k, 2^1(0, 5; c; z) = r(a)r(f)) TT Gauss hypergeometric 

function. For k = n, the partition function C~^ simplifies to 


c-^ = 


n!Pf [7r2crfecr^(crfc - aj)/{ak + f^i)] , 

[7rVfeCrj-(crfe - cr^-)/K + ^ 3 )] 


i!Pf 




n even, 

n odd. 


= n) TT^ 


n 

2=1 j>k 


(7j (7 iz 

O’ 7 + O’Tc 


(90) 


(91) 


The gap probability i?(r, s) is given by Eq. (1^ with 


X3,k{r.s)= J dX J d/r|^[(l + a-2A2)-«+-i(l + a,-V)-"+"-'] 

00 ) ( —oo,r)U(s,oo) 

[ dfx [(1 + a-2A2)-"+"-i(l + (Tfc - (1 + + fTfc 2 a2) 

J /r + A 

(92) 


( —oo,r)U{s,oo) ( —oo,r)U(s,oo) 

5 " 

( —cx3,r)U(s,oo) ( —oo,r)u(s,oo) 


and 


Xj,n+i{r,s) =-Xn+i,j{r,s) = {1-Sj^n+i) J dXil + a^X^) 
^/tt ajT (k — n + I) 


— 2 \2\—K+n—1 


( —cx3,r)u(s,oo) 


3 


3 5^ 


T./ , IN +r 2 Fi\ -,K-n+l;-] - ^ \ - s 2 F 1 \ k - n + 1; 

T[K-n+l) \2 2 aj J \2 2 cr/ 


(1-d,■„+!). (93) 


— K+n —11 

5 
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Figure 5: Plots for correlated Cauchy-Lorentz (Variant I) ensemble with n = 3 ,k = 3.4, and S = (1/7, 3/8, 2/9). 
(a) SF of the smallest eigenvalue (b) PDF of the smallest eigenvalue (c) CDF of the largest eigenvalue (d) PDF 
of the largest eigenvalue. 


Similarly, Xj,k{r, s) as given below yields the double gap probability using Eq. (15^ : 

X,Ar,s)= [(l + a-^A^)-«+"-i(l + a,-V)-'=+"-'] 

= i / dA / dfi [(1 + + (Tfc - (1 + (T-V")“"+”"^(l + CTfc "a2)-"+"-i] , 

^ Jr Jr + A 


(94) 


Xj,n+l{r,s) = -Xn+l,j{r,s) = {1 - Sj^n+i) j dA (1 + cr^- ^A^) 


2 \ 2'\—K+n—1 


1 


3 


1 


3 


I + I I - n + 1; 


(1 ^i,n+l). 


(95) 


The expressions for density of the extreme eigenvalues follow from Eqs. ((36l) and (l38l) . Similar to the correlated 
Gauss-Wigner case, here also we have x(—oo, —x) = xix, oo) which gives E(—oo, —x) = E{x, oo) and ps{—x) = 
Pl{x). 

In Table [S] we compare the gap probabilities obtained using analytical results and numerical simulations. 
Fig. [5] shows the plots of distributions and densities for the extreme eigenvalues. 


19 





































5.2. Uncorrelated case 


The uncorrelated Cauchy-Lorentz (variant I) is obtained from Eq. 
matrices H, 

P(H)cx|1„ + H2|-«, 


for S = We have for hermitian 


(96) 


with K > n — lj2. Interestingly, for n = n the matrices H may be generated using n-dimensional unitary 
matrices U 


H = 


1 /1„-U 


1. 


u 


(97) 


the measure |1„ + being then equivalent to the Haar measure d^(U) Here l = 'J—l is 

the imaginary unit. For arbitrary k, Eq. (1^71) can be used provided U are taken from the measure ||lln, + 
U||^^'^~")d/i(U), where || - || stands for the absolute value of determinant. 

Using Eq. (IMl), because of unitary-invariance, the joint probability density of eigenvalues can be immediately 
written down as 


PiW) = cl[ 


11(1 +A?)" 


A^({A}). 


(98) 


The normalization factor is provided by the result 


C-^ =nl 


1 + (-1)^+'" g p' + fc - 1 ^ i + ^ -1 


_ 2^n'^ —2Kn-\-n^n 


n 

i=i 


2 ’ 2 
r(j + l)r(j + 2 k- 2 n) 

{j + K — n) 




(99) 


B(a, 6) = /p — t)^ = r(a)r(6)/r(a + b) being the Beta function. 

The uncorrelated Cauchy-Lorentz ensemble described by Eqs. (EH) or 
in its eigenvalue spectra. It has been investigated in the context of modeling financial correlations 


is known to exhibit Levy tails 

mMI- In 


the problem of quantum chaotic scattering, Cauchy-Lorentz density has played a crucial role in demonstrating 
equivalence between the Heidelberg (Hamiltonian) approach and the Mexico (Scattering matrix) approach, the 


latter being dictated by the Poisson kernel [90j. Gap probabilities for the uncorrelated case has been studied 
in [l^ and, for finite n case, presented in terms of second-order second-degree ordinary differential equations 
which are related to certain Painleve-VI transcendents. In 0 the density of the largest eigenvalue has been 
calculated for large n and its large deviations have been examined. We provide below explicit finite n results 
as determinantal expressions. 

The expression for gap probability is given by Eq. (1291) with 
(i+(-i)^+U 




+ i(-l)J+feti+fe 


-2k-1 


B(-;I., 1 - «:) - B(- 1 1 - «) 


Xj,k{r,s) = < 


1. ,j-\-k — 2K— 1 
2^ 

(i+(-ip+U 

2 


(_ip+fcB(-A, K _ 1±|^, 1-k) + B(-- 


2 

.i+k-1 


1 — k) 


B(^±|^,K-i±|^) 


.7+^—1 ' 
2 


_ ^ 2k —1 

2^ 


B(—K — ‘AA-A, 1 — k)— B(—■^,K — 


.7+fc-l 


, 1 - k) 


r < s < 0, 
r < 0 < s, 

0 < r < s. 


( 100 ) 
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Table 6: Gap probabilities: Comparison between analytical and simulation results for uncorrelated Cauchy- 
Lorentz (variant I) ensemble. 


n 

K 

r 

s 

E{r,s) 

E{r, s) 

Analytical 

Simulation 

Analytical 

Simulation 

2 

3 

-5 

-0.3 

0.2907 

0.2904 

0.0307 

0.0309 

2 

1.9 

0.1 

4 

0.3159 

0.3112 

0.0555 

0.0551 

3 

3.5 

-1 

0.6 

0.0303 

0.0300 

0.0238 

0.0239 

4 

5.8 

-0.65 

0.65 

0.0025 

0.0028 

0.0036 

0.0038 

5 

6 

-7 

4.7 

0.0000 

0.0000 

0.9194 

0.9197 



Figure 6: Plots for uncorrelated Cauchy-Lorentz (variant I) ensemble with n = 3, k = 3. (a) SF of the 
smallest eigenvalue (b) PDF of the smallest eigenvalue (c) CDF of the largest eigenvalue (d) PDF of the largest 
eigenvalue. 


Here B(z, a, b) = dt ^(1 ~ tY ^ is the incomplete Beta function. We note that when j + k is odd, all the 
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three expressions above become identical. Similarly, for E{r, s) given by Eq. (1321) we have 






.i+k-l 


,1- k)- B(-^,k - 


.i+k-l 


, 1 - At) 


Xj.fc(r,s) = < 


2 ’ ' 


i+fc-1 


_ ^ .j-\-k—2K — l 


(_iy+feB(- J,, At - 1 - At) + B(-^, At - 1 - At) 


1 .j-\-k — 2K—l 

t s'- 


B(— At — 1 — At)— B(— -ijAt — ■^"''2 ^ , 1 — At) 


r < s < 0, 


r < 0 < s, 

0 < r < s. 

( 101 ) 


For r = —ooors = oo these reduce to 
X3,k{-00,x) = Xj,k{x,OG) 

B(i±|^, At - i±|^) - i(-l)J + fcAJ+fc-2«-lB(- 1 , At - i±|^, 1 - At), X < 0, 


^ - ^±1^, 1 - At), 


( 102 ) 


a; > 0. 


and 


(103) 


Xj,k{x,oo) = Xj,fe(-oo,a;) 

1 (_l)a+fe,a+/c-2«-iB(_^, ^ 1 _ ^ < 0, 

+ B(i±|^, At - ^±1^) - iAJ+'=-2«-lB(-^, At - i+l^, 1 - At), X > 0. 

We remark that the kernels Xj,k{r, s) and Xj,k{x, s) can also be represented in terms of Gauss hypergeometric 
function 2E1, however we have opted to present the results in terms of Beta functions because of their simpler 
nature. 

The expression for density of the smallest eigenvalue is obtained as Eq. (1201) with 

2,i+fc-2^2^ _l_ a;2)-K^ j — 

Xj,k{-oo,x), 

Likewise the expressions for the density of the largest eigenvalue is given by Eq. (1^^ with 

^i+As-2/^_|_ 2 \-k j=i 

= ■ 




(104) 


(105) 


Xj,k{x,oo), j^i. 

Once again we have E(—oo,—x) = E{x,oo) and ps{—x) = pl{x), the reason being similar to that in the 
uncorrelated Gauss-Wigner case. 

Table[n]collates the results for gap probabilities for various choice of parameters. Fig.[S]shows the comparison 
between the analytical and simulation-generated plots for extreme eigenvalue distribution and density functions. 


6. Cauchy-Lorentz ensemble: Variant II 

6.1. Correlated case 

Consider n x ua dimensional complex matrices A from the density 

iPA(A) oc |l„-bS-iAAt|-", (106) 
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with UA > n, K > UA + n — 1, and S = diag (tri,tT„), CTj > 0. We are interested in the eigenvalues of 
the positive-definite Hermitian matrices H = AA^, or equivalently the singular values of A. The probability 
density function satisfied by H is 

iP(H) (X (107) 

We many generalize this density and replace ua — n by a real a > —1, so that 

iP(H) (X |H|“|l„ + S-iHr«, (108) 


with K> 2n + a — 1. We refer to this as the correlated Cauchy-Lorentz (variant II) ensemble. 

This ensemble admits construction via ratio of two n-dimensional Laguerre-Wishart matrices, namely 


H = 


Ws' 


where W a and W b are respectively from the densities 


Pa(Wa) cx IW^re- 


Pb(Wb)oc 


Wb 


(109) 


( 110 ) 


such that K = a + /3 + 2n, and S We have assumed here /3 > —1, similar to a. A proof has been 

provided in [Appendix A[ The above construction is particularly useful for generating matrices if a (> — 1) and 
K (> 2n + a — 1) happen to be integers. This is because in this case the n-dimensional Wishart matrices 
and Wb with respective degrees of freedom n^, ns can be easily generated with the aid of Eq. (IMl) . and then 
a = Ua — n, K = Ua + ns- We also note that the matrices b and share 

the same set of eigenvalues as they correspond to the identical generalized eigenvalue problem. Therefore, as 
far as eigenvalue statistics is concerned, we may use the notation in Eq. (I109D without any ambiguity. It is also 
clear from Eq. (II09I) that we are dealing with a multivariate generalization of beta distribution of the second 
kind. 

The joint density of eigenvalues is obtained on employing the unitary group integral result given in as 


P({A}) = C'A({A}) nAM(l + a-'Afc) 


— K;+n—1 1 




( 111 ) 




with 0 < Xj < oo. Clearly we are dealing with type I ensemble as in Eq. dini) with w{X) = A“, /j (Afc) = Af 
and gj{\k) = (1 + 0-“^Afe)“'^+"“^. The partition function is given by 


C-^ = n\ 


B(j + a, K-a-n+l-j) 




= n\ A({(t}) ]q cr“+^ B(j + a,K - a-n+l- j). 


( 112 ) 


i=i 


The above variant of Cauchy-Lorentz model has been used in 


9l| to work out the eigenvalue statistics of 


correlated Jacobi ensemble; see section FTTl Moreover, the JPDE in Eq. (Illip has already been derived therein. 
We explore the behavior of its extreme eigenvalues below. 
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Table 7: Gap probabilities: Comparison between analytical and simulation results for correlated Cauchy- 
Lorentz (variant II) ensemble. The a values are from (cti, 0-5) = (1/2,2/3,7/6,3/10,3). The uajTIb values 
are indicated whenever a straighforward matrix construction using Eq. (11091) is possible. 


n 

a 

K 

UA 

riB 

r 

s 

E{r,s) 

E{r, s) 

Analytical 

Simulation 

Analytical 

Simulation 

2 

1 

5 

3 

2 

0.4 

6 

0.2505 

0.2517 

0.1572 

0.1563 

3 

-0.6 

5.9 

- 

- 

0.23 

5 

0.0844 

0.0822 

0.0032 

0.0032 

3 

0 

7 

3 

4 

0.2 

3 

0.0501 

0.0503 

0.0130 

0.0131 

4 

1 

9.5 

- 

- 

0.2 

5.3 

0.0015 

0.0014 

0.0048 

0.0043 

5 

2 

13 

7 

6 

3 

5.5 

0.5605 

0.5608 

0.0000 

0.0000 


The expression for gap probability E{r, s) is as in Eq. (1141) with 

X 3 ,k{r,s) = [(-1)“^““B (-cr^V,j +a,n- k) - (-l)-^+“+”“'" B (-<7^ k - a - n + 1 - j,n - k)] . 

(113) 

While the double gap probability E{r,s) is obtained from Eq. (ITSl) using 

XjAes) = (-l)"'^"“cr^’^“ [B(-cr^^s,j + a,n-K) - B j + a,n - k)] . (114) 

For r = 0 or s = 00 the above two equations simplify to 


Xi,fe(0,a:) = Xj,fc(a;, 00) =-(-1)-^+“+" B (-crfc s \ k - a - n + 1 - j, n - k) , (115) 


Xi,fc(a:,oo) = Xj.fe(0,a;) = (-1) ^ “cr^’^“B (-cr^, V, j + a, n - k) . 
The expression for density of the smallest eigenvalue is obtained as Eq. (1201) with 


= 




(116) 


(117) 


J = b 

Xi.fc(0,a;), j^i- 

Parallely, the expressions for density of the largest eigenvalue is given by Eq. (1^^ with 

X'’+“-i(l + cr/'ia;)-'^+"-i, j = i, 

^Xj,kix,oo), j^i. 

In Table[7]we present the gap probability values for several combinations of parameters. Fig. [3 exhibits the 
behavior of extreme eigenvalues. 


(118) 


6.2. Uncorrelated case 

We now examine the uncorrelated version (S = !„) of Eq. p08D . namely the density 


iP(H)cx |Hni„ + H|-% 


(119) 
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Figure 7: Plots for correlated Cauchy-Lorentz (variant II) ensemble using n = 4, riyi = 4, = 5 which gives 

TO = 4, a = 9. Also, 'Ea = diag (33, 21,15,43), = diag (2,1, 3,4) which gives S = diag (33/2, 21, 5,43/4). 

(a) SF of the smallest eigenvalue (b) PDF of the smallest eigenvalue (c) CDF of the largest eigenvalue (d) PDF 
of the largest eigenvalue. 


involving positive-definite Hermitian matrices H with k > 2n + a — 1. In accordance with the correlated case, 
this ensemble is related to 

iP^(A)oc |l„ + AAt| —, (120) 

with A being a rectangular complex matrix. The form (11201) of uncorrelated Cauchy-Lorentz ensemble has 


been explored in connection with the projection formula in supersymmetry 


92|. 


Similar to the correlated case, this ensemble can be realized using the ratio of two Wishart matrices, i.e., 
H = Wa/Wb where now and are from the distributions 


Pa(W^) cx |W^re-‘’''^A Pb(Wb) cx 


( 121 ) 


such that K = a + /3 + 2n. As in the correlated case, for integer a (> —1) and k (> 2n -|- a — 1) the above 
construction is particularly useful. The matrix model H = W a /W b follows in a special limit (a = 6 —>■ oo) of 
a more general case aW a/ (In + b'W b)\ a,b> 0 considered in [66j. 

The joint density of eigenvalues follows readily because of unitarily-invariant nature of Eq. (11191) . We have 


P({A}) = CA2({A})^A“(l + A,)-^ 
i=i 


(122) 
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Table 8: Gap probabilities: Comparison between analytical and simulation results for uncorrelated Cauchy- 
Lorentz (variant II) ensemble. The n^, ub values are indicated whenever an easy matrix construction is possible 
using ratio of two Wishart matrices. 


n 

a 

K 

UA 

UB 

r 

s 

E(r, s) 

E{r,s) 

Analytical 

Simulation 

Analytical 

Simulation 

2 

0 

5 

2 

3 

0.9 

5.6 

0.3169 

0.3134 

0.0122 

0.0120 

3 

-0.8 

6.7 

- 

- 

0.3 

11.3 

0.0157 

0.0159 

0.0011 

0.0010 

3 

1 

9 

4 

5 

0.6 

24 

0.0028 

0.0025 

0.0096 

0.0092 

4 

1 

11 

5 

6 

3.4 

CX) 

0.2118 

0.2111 

0.0000 

0.0000 

5 

2 

13.7 

- 

- 

0 

3.5 

0.0000 

0.0000 

0.0155 

0.0155 



Figure 8: Plots for uncorrelated Cauchy-Lorentz (variant II) ensemble using n = 2, ha = 7, ub = 2 which gives 
q: = 5, K = 9. (a) SF of the smallest eigenvalue (b) PDF of the smallest eigenvalue (c) CDF of the largest 
eigenvalue (d) PDF of the largest eigenvalue. 


The normalization factor is obtained using 


C ^ = n\ |B(j + k + a-l,K-a-j-k + l)|j,fe=i,...,„ 

, IT r(j)r(j + a)r(K -a-n- j + l) 

= "'0 --• 


(123) 


26 



































































































































The gap probabilities are given by Eqs. da cinci (jlSp with 


s) =-(-1) ^ “ B (-r, a + j + A: - 1,1 - k) - '"B(-s \ k - a - j - fc + 1,1 - k) , 

(124) 

and, 


XjAr,s) = (- 1 ) 


—j — k — cx. 


B (—r, a + j + A: — 1, 1 — At) — B(—s, a + j + A: — 1, 1 — k) 
These simplify to the following for r = 0 or s = oo: 

Xj,fe(0,x) = XjAx,oo) = -(-l)^+''+““'^B (-a::"\K - a - j - A: + 1,1 - At) , 


(125) 


(126) 


Xj,k{x,oo) = Xj,kiO,x) = -{-!) ^ “B(-a;,a + j + A:- 1 , 1 -At). 


(127) 


The expression for density of the smallest eigenvalue is obtained as Eq. ((^(1 with 


4>%{x) = 


^a+j+k 2 _|_ 2 ;) K.^ j = i, 


Xi,fc(0,a;), j^i. 

while, the result for density of the largest eigenvalue is given by Eq. (l22l) with 

,,b) 


AA ^) = 


^a+j+k 2 j = i, 

Xi,fc(a;,oo), j^i. 


(128) 


(129) 


In Table [8] we collect the gap probability values for various choice of parameters. Eig. [8]shows the plots 
of distributions and densities obtained from analytical expressions (solid lines) as well as from Monte-Carlo 
simulations (histograms). 


7. Jacobi-MAN OVA ensemble 

7.1. Correlated case 

Consider the matrix probability density 

iP(H) oc |Hr \tn - nf |1„ + (S”1 - 1„)H|-", (130) 

with H being n x n-dimensional matrices, such that 0 < H < and S = diag (cri,..., (t„), Uj > 0. By the 
notation H > K for matrices, we mean that H — K is non-negative definite. Also a, [3 > —1 and k is any real 
parameter. Eq. (11301) then dehnes the correlated Jacobi-MANOVA ensemble of random matrices. The term 
MANOVA derives from the area of multivariate statistics and stands for Multivariate Analysis Of VAriance, 
while Jacobi has to do with the occurrence of Jacobi weight function in the JPDF of eigenvalues; see below. 
We note that the eigenvalues for matrix model (11301) varies from 0 to 1. In contrast, if we consider the density 

nu) cx |i„ -b Hr \in - Hr i(s-i +1„) -b (s-1 - i„)H|-r (m) 


27 








then the corresponding eigenvalues range from —1 to 1. The two models, (11301) and (I131|) . are trivially related 
by a shift and a scaling transformations, and so are the corresponding eigenvalues. 

For P = K — a — 2n > —1 the correlated Jacobi-MANOVA ensemble has a direct relationship with the 
Cauchy-Lorentz (variant II) ensemble, Eq. (11081) . For if we consider matrices H drawn from Eq. (|108l) . then 
the matrices 


H = ^ (132) 

In+H 

are distributed according to the density (11301) . Equivalently, if we consider matrices the matrices H drawn 
from Eq. (I130L then the matrices 

H 


«=wi' 

are distributed according to the density (11081) . As a consequence of this correspondence, we may construct 
the Jacobi-MANOVA matrices as follows. Consider n-dimensional Laguerre-Wishart matrices and 
drawn from the densities as in Eq. (IllOl) . then 


H = 


W, 


(134) 


Wa+Wb 

constitutes the ensemble described by Eq. (I130L with a,/3 as in Eq. (IllOp and k = a + P + 2n. Also, 
s-i = In [Appendix B| we provide a proof. Again, for (permitted) integer values of a, p, k the 

above construction gives a very accessible way of generating Jacobi-MANOVA matrices, as the n-dimensional 
Wishart matrices W^, with respective degrees of freedom n^, us can be easily constructed using Eq. (|M1). 
In this case the parameters are given by a = ua — n, P = ub — n and k = ua + tib- From Eq. (11341) it is also 
clear that we are dealing with a multivariate generalization of beta distribution of the first kind. Correlated 
Jacobi-MANOVA ensemble has been introduced very recently and solyed exactly to obtain the joint eigenyalue 


density and correlation function of arbitrary order 


A unitary group integral result from 
density (jl30p as 


91|. 


leads to the joint probability density of eigenyalues for matrix 


P(Ai,..., A„) = CA({A}) n Ar(l - [1 + (^7' - l)Afe] 


1^1 


] —K-\-n— 1 




(135) 


Note that for k — n-l-1 = 0, one encounters 0/0 form in the JPDF of eigenyalues as well as in the expressions for 
gap probabilities and extreme eigenyalue densities, therefore a limiting procedure must be inyoked. Howeyer, 
this case (along with S = In) corresponds to the uncorrelated Jacobi-MANOVA ensemble for which much 
simpler results exist, as giyen in the next section. Another case when a limiting procedure has to be followed 
occurs when k — n-|-l = —1 with n > 2. 

As a consequence of the relationship with Cauchy Lorentz (yariant II), when P = k — a — 2n the results for 
correlated Jacobi-MANOVA can be written down easily by a simple yariable transformation. This also includes 
the case when construction of Jacobi-MANOVA matrices is possible using two Laguerre-Wishart matrices, as 
described aboye. Therefore, this scenario is of special interest and we explicitly deal with this case along with 
the general k case. 
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7.1.1. Arbitrary k, and a, j3 > —1 

The partition function is obtained as 

= n! |B(j + a,/3 + 1) 2 -Fi(j + a, K - n + 1; j + a + /3 + 1; 1 - (136) 

where 2 ^ 1 ( 0 , b; c; z) is the Gauss hypergeometric function, as already mentioned. 

The gap probability is given by Eq. (HI with 


Xj,k{r,s) = B{j + a,/3 + 1) 2 Ei(j + a, K - n + l;j + a + /3 + 1; 1 - cr^ 




+ — ^1 (j + a; -/3, K - n + 1; j + a + 1; r, (1 - ^)r) 

J +a 




■ Fi [j + a; -/3, k - n + 1; j + a + 1; s, (1 - cr^ ^)s) 


r(a+m+ra)r(bi+m)r(fe 2 +n)X y j-gpj-gggj^^g ^j^g Appell hyperge- 


j + a 

where Fi(o; 61 , 62 ; c; x, y) = r(a)r^Sr(b 2 ) Sm,n=o- m\n\r{c+^+u) 

ometric function of two variables. Similarly, we have 

gj+a 

Xj,k{r,s) = — — El {j + a;-^,K-n + 1; j + a + l;s, (1 - cr^^)s) 
j + a 


(137) 




J + a 

For r = 0 or s = 1, we have the following simplifications: 


Fi {j + a; -P,K- n + 1 ; j + a + l;r, (1 - cTfc ^)r) . 


(138) 


Xj,k(.0, x) = Xl.k{x, 1) = B(j + a,j3 + l) 2 F 1 (j + 
x3+°‘ 

- --— ^1 (j + a; -/3, K - 

3 + a 


Of, K — n + 1; j + a + /3 + 1; 1 — (T^ 
n + 1 ; j + a + 1 ; a;, (1 - iyk^)x) , 


(139) 




X3,k{x,l) =Xj,k{^iX) = -—Fi {j + a-,-l3,K-n + l\j + a + l]x,{l-(jf.^)x) 
J “r Oi. 


The density of smallest eigenvalue is given by Eq. (uni) with 

a:l+“-i(l — a;)^(l + (cr^^ — 1 ) 2 :)“''+"'“^, j = i. 




3 7^ b 


^Xi,fc(0.a:), 

and that of the largest eigenvalue by Eq. (1^^ with 

a;l+“-i(l — a;)^(l + (cr^^ — l)a;)“'^+"“^, j = i, 




X],k{x, 1), 


3 ^ i- 


7.1.2. Arbitrary a > —1, and fd = n — a — 2n > —1 


In view of the result (|132l) we consider the transformation 911 


Afe = —-, /c = l,...,n. 

1 + Mfc 


(140) 


(141) 


(142) 


(143) 
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Table 9: Gap probabilities: Comparison between analytical and simulation results for correlated Jacobi- 
MANOVA ensemble. The a values are from (cti, 0 - 5 ) = (3/2, 5/6, 7/4, 2/5,4). The ha, ns values are in¬ 
dicated whenever a straightforward matrix construction using Eq. (I134p follows. 


n 

a 


n 

riA 

riB 

r 

s 

E{r, s) 

E{r, s) 

Analytical 

Simulation 

Analytical 

Simulation 

2 

1 

2 

7 

3 

4 

0.2 

0.6 

0.3128 

0.3133 

0.1071 

0.1065 

2 

- 0.1 

- 0.3 

2 

- 

- 

0 

0.5 

0.0891 

0.0883 

0.0352 

0.0356 

3 

1.3 

1.4 

1.5 

- 

- 

0.5 

1 

0.0029 

0.0026 

0.0034 

0.0035 

4 

1 

2 

11 

5 

6 

0.1 

0.7 

0.0013 

0.0013 

0.0058 

0.0054 

5 

2.2 

4.5 

0 

- 

- 

0.15 

0.9 

0.0000 

0.0000 

0.0702 

0.0704 


The joint density in the new variables is 

n 

P(^i ,...,= CAiif,}) (1 + a- 


Z=1 




Here the normalization factor C is same as the C in Eq. (11121) . Also, C of Eq. (I136p reduces to C when 
f3 = K — a — 2n. We use Eqs. (11121) — (jllSp to write the expressions for gap probabilities and extreme eigenvalue 
densities in the /i variables, and then transform back to the A variables. We wind up with the results as given 
below. 

The gap probability is given by Eq. (HI with 


Xj.fe(As) = cr; 


J+a 


(_1)-1-«B 


(Tk(r - 1) 


, /-|-Qf,n— k) —(—1)1+“+" B ( - 11 

(144) 


While, the double gap probability is obtained using Eq. (fT51) with 


X,,fc(r,s) = (-l)-^-“ai+“ 
For r = 0 or s = 00 these reduce to 


B 


tTfc(s - 1) 


, j + a,n — K \ — B 


o-fe(r - 1) 


,j + a,n — K 


(145) 


Xj.fe(0,a:) = Xj,fc(x,oo) = -(-1)^+“+" B fZhil —ll^^- a - n-bl-j,n-«:), (146) 


(147) 


= Xi,fc(0,a;) = (-1) ^ “cr^+“B^- 


-,j +a,n- k] . 


^ak{r - 1)' 

The probability density function for the smallest eigenvalue is obtained from Eq. dsni) after incorporating 
Jacobian of transformation as 


psix) = n\C '^\<t>%{x)\j,k=i,...,n, 

^ ' i—1 


(148) 
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Figure 9: Plots for correlated Jacobi-MANOVA ensemble. (a),(e) SF of the smallest eigenvalue; (b),(f) PDF 
of the smallest eigenvalue; (c),(g) CDF of the largest eigenvalue; (d),(h) PDF of the largest eigenvalue. For 
(a)-(d), the parameter values used are n = 3,a = —0.2,/3 = 1,k = 2.5. Also, S = (4,7/5,10/3). This set of 
values does not admit a straightforward matrix construction using Eq. (11341) , and therefore results from section 
7.1.1 have been used for the analytical plots. For (e)-(h) the parameter set is such that matrix construction 
using Eq. (jl34l) is easily possible and analytical results from section 7.1.2 have been used. The parameter 
values are n = 4,nA = G,nB = 5, and = (5/2,1/11,2/7,1/4), = (7/6,1/5,4/3,1/3). These lead to 

a = 2,/3 = 1,K = 11, and S = (15/7,5/11,3/14,3/4). 
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with 





1 — X J 





_^) 

(Tfe(l - X) J 


Xj.fe(0,a:), 


The expressions for density of the largest eigenvalue is given by 


, J = *, 

j *• 


Pi_{x) =n\C 


1 

{i-xY 


2=1 


(149) 


(150) 


with 





x 

(Jk{l - x) 




j = i, 


(151) 


[xi,fc(a;, 1), j^i- 

In Table [3] we compare the gap probability values obtained from analytical formulae and from numerical 
simulations. Fig. IH] shows the plots of distribution functions and probability density functions. For (a)- 
(d), we have used the parameter values n = 3, a = —0.2,/3 = 1 ,k = 2.5. Also, the correlation matrix 
taken is S = (4,7/5,10/3). This set of values does not admit matrix construction using Eq. (11341) in a 
straightforward manner, and therefore we use the results from section 7.1.1 for the analytical plots. For (e)-(h) 
the parameter set is such that matrix construction using Eq. (11341) is trivially possible and we use the analytical 
results from section 7.1.2. The parameter values are n = 4, ua = 6,nB = 5, and Sa = (5/2,1/11,2/7,1/4), 
Ss = (7/6,1/5,4/3,1/3). These lead to a = 2,/3 = 1, k = 11, and S = (15/7,5/11,3/14,3/4). We find perfect 
agreement between the analytic predictions and the simulation results. 

We note that E(0,x), E(x, 1), ps(x) and pl(x) for a given set of n,a,/3, S are, respectively, same as 
^(l — X, 1), £'(0,1 — x), pl(1 — x) and ps(l ~ x) for n,/3,a, This correspondence has to do with the 

structure of the matrix density (11301) . 


1.2. Uncorrelated case 

The uncorrelated Jacobi-MANOVA case follows from Eq. (11301) for S = 1„. We have 

£(H)oc |Hni„-H|^ (152) 

where a, (3 > —1. The eigenvalues of this model are from 0 to 1. On the other hand if we work with the model 

£(H)cx|l„ + Hni„-H|^ (153) 

then the eigenvalues he in [—1,1]. Again, the two matrix models are related by a simple linear transformation 
and so are the eigenvalues. As in the correlated case, the matrices from (11521) can be generated using n- 
dimensional Laguerre-Wishart matrices as H = Wa/(Wa + Wb) where Wa, Wb are from (I12ip . 

The joint eigenvalue density is obtained as 

n 

P(Ai,..., A^) = CA2({A}) n A“(l - A,)^. (154) 

i=i 
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Table 10: Gap probabilities: Comparison between analytical and simulation results for uncorrelated Jacobi- 
MANOVA. The nA,nB values are indicated whenever an easy matrix construction is possible using H = 
W^/(W^+Wb). 


n 

a 

/3 

riA 

riB 

r 

s 

E{r, s) 

E{r,s) 

Analytical 

Simulation 

Analytical 

Simulation 

2 

3 

1 

5 

3 

0.3 

0.7 

0.1977 

0.1970 

0.0850 

0.0855 

3 

-0.4 

3.5 

- 

- 

0.2 

0.8 

0.0179 

0.0178 

0.0029 

0.0028 

3 

0 

4 

3 

7 

0.05 

0.4 

0.0555 

0.0538 

0.0146 

0.0141 

4 

3.2 

1.8 

- 

- 

0 

0.37 

0.0233 

0.0232 

0.0000 

0.0000 

5 

2 

9 

14 

7 

0.5 

1 

0.0000 

0.0000 

0.0135 

0.0130 


This JPDF is trivially related to the standard Jacobi ensemble of random matrices 
transformation Xj —^ (1 + Xj)/2. The normalization factor is obtained using 


E 


63|, |77| via the linear 


C — n! |B(j + A: + a — 1,/3 + 1) 


n 

i=i 


r(j + i)r(j + a)r(j + /3) 


r(j + a + /3 + n) 

Jacobi-MANOVA ensemble has an importa nt ro le in multivariate statistics 


(155) 


tions in quantum conductance problem [& 


781, 


79( 1 ■ Moreover, it fin d applica 


and multiple channel optical fiber communication 


103 


The extreme eigenvalue statistics for uncorrelated Jacobi-MANOVA has been investigated for large n in 


28 


104|. 


33|. 


Finite n case has been considered in [29j. However, even for the complex case the results have been provided 
in terms of hypergeometric function of matrix argument. We provide below expressions which are in terms of 
standard Beta functions, and are much easier to implement. 

The expression for gap probability for the eigenvalues described by Eq. (11541) is given by Eq. (|29|) with 


Xj,k{r, s) = B(j-|-fc-|-a—l,/3-|-l)-|-B(r;j-|-/c-|-Q!—l,/3-|-l) — B(s; j + k + a — l,/3-|-l). (156) 


Similarly, 


Xi.fc(As) = B(s;j + k + a-1,13+1)- B(r;j + k + a-1,13 + 1) 


(157) 


leads to the expression for double gap probability using Eq. These simplify to the following for r = 0 or 

s = 1: 

Xj,ki0,x) = X 3 ,k{x, 1) =B{j + k + a -l,j3+ 1) - + k + a - 1, (3 + 1), (158) 


Xj,/c(a:, 1) = Xj,k{'3,x) = B(x;j + k + a - 1,13 + 1). 


(159) 
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Figure 10: Plots for uncorrelated Jacobi-MANOVA ensemble with n = 5,nA = 7, ns = 8. (a) SF of the 
smallest eigenvalue, (b) PDF of the smallest eigenvalue, (d) CDF of the largest eigenvalue, (d) PDF of the 
largest eigenvalue. 


The probability density function for the smallest eigenvalue is given by Eq. (120)1 with 




1 +j+fe j = i, 


(160) 


Xj.fc(0,a:), j^i. 

Likewise, the expression for the probability density of the largest eigenvalue is obtained as Eq. 

n 

pi^{x) =n\c'^ IV’j*i(a;)|j,/c=i.....n, 


with 


(161) 


where 


j ^a+j+k 2^2^ _ 2-^/3^ j = i, 

1), jV*- 


(162) 


In Table (TU] we compare the results for gap probabilities for various choice of parameters. Fig. [TU] shows 
the plots of extreme eigenvalue distributions and densities obtained from analytical expressions as well as from 
Monte-Carlo simulation. 

Similar to the correlated case, E{0,x), E(x, 1), ps{x) and PLix) for a given set of n,a,/3 are, respectively, 
same as ^(l — x, 1), E{0 ,1 — x), pl(1 — x) and ps(l — x) for n, /3, a. 
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8. Bures-Hall ensemble 


8.1. Correlated case 

Bures-Hall measure u sually pertains to a fixed trace scenario, its origin being in the field of quantum 
information 


87 


105 


- lOSj . However, we retain the name in the case of unrestricted trace also. We can write 
down the matrix probability density for positive-definite Hermitian matrices H as 

1 


iP(H) oc |H|“e 




(163) 


A+(H)’ 

where a > —1 is a real parameter, S = diag (oi,..., cr„), Oj > 0, and A-|_(H) = A+({A}), A’s being the 
eigenvalues. We can realize the above ensemble using 

W 


where the W’s are n-dimensional matrices from the complex Wishart distribution P(W) oc |W|'’ exp(— tr S~^W), 
as in Eq. dm]), but with rj > —1/2. While the G’s are n-dimensional Hermitian matrices from correlated 
Cauchy-Lorentz distribution (Variant I), 

iP(G) oc |11„ + SG2|-(’7+"), (165) 

cf. Eq. (1861) . The parameter a of Eq. (I163p is then given by a = 77 — 1/2. A proof has been outlined 
in [Appendix C[ A special case occurs when the W are n-dimensional Wishart matrices with degree of freedom 
TO, so that f] = m — n. In this case a (> —1) assumes half-integer values, namely a = —1/2,1/2,3/2 etc. 

The joint density of eigenvalues turns out to be 


P({A}) = C. 


A({A}) 








(166) 


a+({a})ia- 

with Aj > 0. This JPDE is of type H since it can be rewritten as the product of a Pfaffian and a determinant, 
as in the correlated Gauss-Wigner case. The partition function is obtained as 

C~^ = n!Pf[hjk]j,k=i,..,N, 

where N is as defined in Eq. (1^51) . Eor even n we have 

hj,k — —^ 27^1 (2o -|- 2, a + 2; 2 q; -|-3;1 — 2 F 1 (20 -I- 2, o -I- 2; 2o -I- 3; 1 — Cj ^o^)] . 

(167) 

and when n is odd, in addition we have 

h,,„+i = = af+^Tia + 1)(1 - (5,.„+i). (168) 

Eor o = —1/2, C~^ simplihes to 

n! Pf - ^)/’ 

[7ryo^^(yb^ - ^)/(^ -F ^)] . ^ 


n even, 


c-^ = 


n\ Pf 


= n! 


n 


2=1 


j>k 


[- 


/ ^j y/Fk 


fc— 


0 


n odd. 


(169) 


(170) 
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Table 11: Gap probabilities: Comparison between analytical and simulation results for correlated Bures-Hall 
ensemble. The a values are from (cti, (T 5 ) = (3/7,8/9,7/11,2/5,3). The m values are indicated when W in 
Eq. (11641) are Wishart matrices with degree of freedom m. 


n 

a 

m 

r 

s 

E{r,s) 

E{r, s) 

Analytical 

Simulation 

Analytical 

Simulation 

2 

- 0.5 

2 

0.5 

5 

0.2263 

0.2268 

0.0229 

0.0216 

2 

3 

- 

2.9 

7 

0.2603 

0.2605 

0.0153 

0.0158 

3 

3.5 

7 

0 

3.2 

0.0011 

0.0011 

0.0134 

0.0138 

4 

0 

- 

1.3 

5.2 

0.0511 

0.0530 

0.0000 

0.000 

5 

2.3 

- 

0.6 

11.7 

0.0000 

0.0000 

0.0361 

0.0368 



Figure 11: Plots for correlated Bures-Hall ensemble with n = 5, a = 16/7 and S = (4/11,3/4,4/3,8/5,9/4). 
(a) SF of the smallest eigenvalue, (b) PDF of the smallest eigenvalue, (d) CDF of the largest eigenvalue, (d) 
PDF of the largest eigenvalue. 
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Gap probability E(r, s) is given by Eqs. (l29l) with 


Xj,k{r,s)= J d\ J 


II A 

(0,r)U(s,oo) (0,r)U(s,oo) 

j iX j 

(0,r)U(s,c>D) (0,'r)U(s,oo) 


and 


Xj,n+l('^;-5) — s') — (1 ^_7,n+l) J dX X e ^ 

(0,r)U(s,oo) 

= cr“+^ [7(0 + 1, a~^r) + r(a + 1, (Jj^s)] (1 - 
The double-gap probability E{r,s) is given by Eq. (15^ with 

Xj,k{'^^s)= [ d\ [ 

Jr Jr li ~T A 

II — X 




-e i '^e >= 


and 


Xj,n+i{r,s) =-Xn+I,j{r,s) = {1- 6j,n+i) / dAA“e ^ 

J r 

= cr“+^ [r(a -h 1, cr“V) - r(a -t 1, cr“^s)] (1 - i5y„+i). 


(171) 


(172) 

(173) 


(174) 


(175) 


In Table fTll we list the gap probability values obtained from analytical results and from numerical simulation. 
Eig. 111! shows behavior of extreme eigenvalues in terms of distributions and densities. 


8.2. Uncorrelated case 

The matrix probability density in this case is 


P(H) oc |H|' 




1 


(176) 


A+(H)’ 

for positive-definite Hermitian matrices H, with a > — 1. Similar to the correlated case we can associate to 
Eq. (I176|) the matrix model 

where W is an n-dimensional uncorrelated Laguerre-Wishart matrix as in Eq. (123 with the replacement 
a —>■ ?7 (> —1/2), and G is n-dimensional Hermitian matrix from uncorrelated Cauchy-Lorentz (variant I) 
density, i.e., from Eq. (I165p with S = The parameter a of Eq. (11761) is then given hy rj — 1/2. Again, if 
the W’s are n-dimensional Wishart matrices with degree of freedom m, then a assumes the half-integer values 
decided by to — n — 1/2. In this case, using the results in we know that H can also be generated using 
W’s and n-dimensional unitary matrices U from the measure ||ll„ -I- U| |^*-"*~”^c?/j,(U) as 
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Table 12: Gap probabilities: Comparison between analytical and simulation results for uncorrelated Bures-Hall 
ensemble. The m values are indicated when W in Eq. (I177p are Wishart matrices with degree of freedom m. 


n 

a 

m 

r 

s 

E{r,s) 

E{r, s) 

Analytical 

Simulation 

Analytical 

Simulation 

2 

-0.5 

2 

0.3 

4 

0.0925 

0.0920 

0.1314 

0.1311 

3 

0 

- 

1 

10.5 

0.0049 

0.0048 

0.0098 

0.0093 

3 

1.7 

- 

2.3 

8.3 

0.0586 

0.0574 

0.0078 

0.0072 

4 

1.5 

6 

1.4 

8.3 

0.0014 

0.0013 

0.0046 

0.0047 

5 

3.5 

9 

8 

12 

0.3019 

0.3015 

0.0000 

0.0000 


For m = n, i.e. a = —1/2, the measure ||1„ + U|simplifies to d^(U), which is the Haar measure 
on the group of n-dimensional unitary matrices. 

The joint density of eigenvalues is obtained as 


where the partition function is given by 

k- j 


c-^ = 


i! Pf 


r! Pf 


^/2 


j + k + 2a 
k-j 


r(j + Q;)r(fc + a) 
r(j + Q;)r(fc + a) 




j + k + 2a 

[-r(fc + a)]fc=i 

r 0 + i)r(j + 2a + i) 


[r(j + a)].^i 




n even, 


odd. 


(180) 


2n2+2o 


n 


r(j + a + 1/2) 


(181) 


fl¬ 


it is 


The above ensemble has been recently shown to be related to the Cauchy two-matrix ensemble in 
worth a mention here that the normalization factor for fixed trace Bures-Hall ensemble had remained elusive 
for long 


105, 

106 

1 , and was finally worked out in 

107 


lated Bures-Hall ensemble has a simple relationship (via Laplace transform) with the normalization factor for 
unrestricted trace uncorrelated Bures-Hall ensemble, the above result readily leads to the former. 

The gap probability is given by Eq. (1^ with 


Xi.fc(cs)= J dX J 


dfjL 


(0,r)u(s,oo) (0,r‘)U(s,cx:)) 

w - / 

(0,r)U(s,c>o) (0,r)U(s,oo) 


dfi 


fji X 




/i -f A 

„-A„- 


(182) 
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Figure 12: Plots for uncorrelated Bures-Hall ensemble with n = 4 and a = —1/2. (a) SF of the smallest 
eigenvalue, (b) PDF of the smallest eigenvalue, (d) CDF of the largest eigenvalue, (d) PDF of the largest 
eigenvalue. 


and 


Xj,n+i(r,s) = -Xn+ij(r,s) = (1 - ^i.n+i) J dA A“+-^ ^ 

(0,r)U(s,oo) 

= bU + a,r) + F(j + a, s)] (1 - (5j>+i). 
The double-gap probability is given by Eq. (l32ll with 
XjAr.s)= f d\ 

Jr Jr /i H” A 

= 11 dX [ - AA) , 

2 J j> J y. ^ “h A 


and 


Xi.n-n(Ds) =-Xn-Hij(r,s) = (1 - / dAA“+^ ^ 

J r 

= [r(j -ba,r) - r(j -ba,s)] (1 - (5j>+i). 


(183) 


(184) 


(185) 


In Table [12] we compare the results for gap probabilities for several parameter choices. Fig. [12] shows the 
plots for distribution functions and the densities for the smallest eigenvalue as well as the largest eigenvalue. 
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9. Summary and Outlook 


In this work we derived exact results for gap probabilities and probability densities of the extreme eigenvalues 
of six matrix ensembles which are of central importance to random matrix theory. We considered both correlated 
and uncorrelated variants of these ensembles. In the process we also proposed generalized versions of some of 
these ensembles and came across some interesting relationships. All the analytical results for gap probabilities 
and probability densities were validated by numerical simulations based on exact diagonalization of matrices 
or Dyson’s log-gas method for generating eigenvalues. 

Our focus here has been on ensembles of complex random matrices. It is known that in the rotationally 
invariant cases the ensembles comprising real a nd quater nion matrices admit joint eigenvalue densities which 


involve product of a determinant and a Pfaffian 


log 


n3- 


The antisymmetric kernels in volved in the Pfaffian 


being respectively sgn(A—^)/2 and in the orthogonal and symplectic ensembles |109l4lll| . Therefore, 

it is possible to write down the results for gap probabilities and densities of extreme eigenvalues for these 
ensembles using the type II ensemble results of section [2l3l As a matter of fact, even for orthogo nal-unitary 


and symplectic-unitary crossovers the joint density of eigenvalues has the type II structure 


63, 


109 


noi, 


therefore results for extreme eigenvalues, in principle, can be worked out. However, in the case of noninvariant 
ensembles involving real or quaternion matrices, unavailability of group integral formulae similar to those in 
the unitary group leads to serious complications in obtaining exact results. In these cases one has to deal 
with hypergeometric functions of matrix arguments, Jack or Zonal polynomials. A more tractable approach 
is based on supersymmetry, using which explicit answers for extreme eigenvalues have been obtained for the 

Moreover, very recently it has been applied to obtain joint density 


correlated Laguerre-Wishart case 
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and correlation functions for eigenvalues of the correlated Jacobi ensemble 


9l|. It will be an interesting but 


challenging task to explore using supersymmetry exact solutions for the real or quaternion counterparts of the 
noninvariant ensembles (other than Laguerre-Wishart) considered in this work, and look for the corresponding 
large-n universal results. 


Appendix A. Derivation of matrix probability density for correlated Cauchy-Lorentz ensemble 

We derive below the matrix probability density for correlated Cauchy-Lorentz ensemble starting from the 
matrix model (11091) . The required density can be obtained as 

P(H)= J d[W^] J d[WB]J(H-(A.l) 

The delta function with matrix argument in the above equation represents the product of delta functions 
with scalar arguments, one for each independent real and imaginary component of the matrix argument. Also 
0(A) represents the matrix theta function, and requires the matrix A to be positive definite (A > 0) or more 
generally nonnegative-deflnite (A > 0) for a non-vanishing result. 
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Employing the transformation W/i —^ W^W^, we get 

P(H)= J d[WA] J d[WB]5(H-WA)|WA|“e-“^(^^''^^+^B')WB|Ws|“+/5+"0(WA)0(WB) 

= |H|“0(H) yd[WB]e“*’'(^^'“+^e')'^®|WBr+'^+”0(WB) (A.2) 

In the second line above, we performed the trivial delta function integral over A. Now B-integral can be easily 
performed to yield 

P(H) oc + s-i|-(«+/3+2’^) 0(H) 

oc |H|“|1„ + 0(H), (A.3) 

which is same as Eq. (11081) . 


Appendix B. Derivation of matrix probability density for correlated Jacobi-MANOVA ensemble 

The density of H as in Eq. (|130l) can be derived using 

P(H)x J d[WA] J d[WB]<5(H-WA(WA + WB)-i)|WAre-*''^^''^^|WB|^e-‘’'^B^WB 

x0(Wa)0(Wb). (B.l) 


Employing the transformation Wb —^ WbW^ leads to 

P(H)x J d[WA] J d[WB]5(H-(!„ + WB)-^)|WAr+'^+”e-*"^A"w.4|w^|/3g-trSeiwBW^ 


X 0(W^)0(Wb). 


(B.2) 


The W^-integral gives 

P(H)x y’d[WB]<5(H-(l„ + WB)”^)|WB|^|S^^ + SB^WBr(“+^+"”^0(WB). (B.3) 

Taking (1„ + Wb)~^ = Y yields 

P(H) X J d[Y] ,5(H - Y)|Y|-2 -|Y-i - - ln)r(“+'^+""^0(Y)0(Y-i - !„), (B.4) 

and then we have the trivial delta-function integral giving 

P(H) X |H|-2"|H-i - + Sb'(H-i - 11„)|-(“+'5+2-)0(H)0(H-i - 1„). (B.5) 

Note that H~^ — 1„ > 0 implies H < l^,. Now, carrying out some readjustments we obtain 

P(H) X |Hni„ - + (S^i - Sb1)H|-(“+'5+2")0(H)0(11„ - H) 

X |Hni„ - H|'5|11„ + (S-i - 1„)H|-(“+'5+2'*)0(H)0(11„ - H), (B.6) 

where = SbS)(^. 


41 






Appendix C. Derivation of matrix probability density for correlated Bures-Hall ensemble 


We derive here the density for correlated Bures-Hall ensemble. Introducing the Fourier representation for 
delta function (5(H — W(ll„ -b using a matrix K, we have 

P(H)oc J d[K] j d[G] j -t-SG^r(’'+")e““'^^''^|W|’'0(W). 

(C.l) 

Integral over W leads to 


P(H)oc J d[K] J d[G]e**’'^“|l„ + SG2r(’’+”) |S-^-bi(l„-bSi/2G^S^/2)-iK|-(''+’"). 

Considering G —^ and employing some readjustments we obtain 

P(H)oc J d[K\ J d[G]e**''^“|S~i-f G2|-(''+”)|l„-bz(S~i-bG^)-iK|-(’7+”^ 

We now use the transformation K ^ (S-1 + G2) K and get 

P(H)oc J d[K\ J d[G]e**''^“(^“'+°")|S”i-bG2|-’’|l„-biK|-(’7+"). 


(C.2) 


(C.3) 


(C.4) 


Ingham-Siegal-Fyodorov formula 112j| can be used to solve the K-integral, and on simplification we have 


P(H)oce-*''^ '“|H|''0(H) J d[G]e-*''®'“. 

The G-integral yields result which depends only on the eigenvalues of H as [8 

/. T n 

1 TT .-1/2 


d[G] 


-trG^H 


A+({A}) 




(C.5) 


(C.6) 


This can be proved using the result from correlated Gauss-Wigner case, with the replacements H —^ G and 
S“^ H in Eq. (I40F We have, with gj representing the eigenvalues of G, 


d[G]. 


-trC'^H 


oc J dgi--- J dgnA'^{{g}) J d/i(U)e 


-trHU+G^U 


OC / - / dgnA'^{{g}) 


1 




A({A})A({32}) 


oc . / dgi--- [ dgn-^^^j^^\e 


(C.7) 


A({A})y J "'‘A+(M)' 

We used above HCIZ formula to perform the unitary group integral. Now, the integrals over gj follow using 
the normalization factor result in Eq. (PI). We have 

1 


d[G]e 


-trG^H 


OC 




j>k ^ 


^(W)/=V 

which gives on little simplification Eq. (jC.6|l . Therefore, the matrix probability density for H may be written 
as in Eq. (jl63ll with the identification a = g — 1/2. 
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